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We propose a non-lattice simulation for studying supersymmetric matrix quantum mechanics in 
a non-perturbative manner. In particular, our method enables us to put M theory on a computer 
based on its matrix formulation proposed by Banks, Fischler, Shenker and Susskind. Here we 
present Monte Carlo results of the same matrix model but in a different parameter region, which 
corresponds to the 't Hooft large-A^ limit at finite temperature. In the strong coupling limit the 
model has a dual description in terms of the DO-brane solution in lOd type II A supergravity. Our 
results provide highly nontrivial evidences for the conjectured duality. In particular, the energy 
(and hence the entropy) of the non-extremal black hole has been reproduced by solving directly 
the strongly coupled dynamics of the DO-brane effective theory. 
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1. Introduction 



Large-A^ gauge theories are playing more and more important roles in theoretical particle 
physics. In particular, they are considered to be useful in formulating superstring/M theories non- 
perturbatively extending the idea of matrix models, which was successful for string theories in 
sub-critical dimensions. For instance, it has been conjectured that critical string/M theories can be 
formulated in terms of matrix models, which can be formally obtained by dimensionally reduc- 
ing U(A'^) super Yang-Mills theory in ten dimensions to D = 0, 1,2 dimensions. The D = 1 case 
corresponds to the Matrix theory which is conjectured to describe M Theory microscopically. 

Another important conjecture, which has been studied intensively over the decade, concerns 
the duality between the strongly coupled \wge.-N gauge theory and the weakly coupled super- 
gravity. The best understood example is the AdS/CFT correspondence, but there are numerous 
extensions to non-conformal field theories as well. In particular, \wg&-N gauge theories in low di- 
mensions have been studied intensively at finite temperature, which revealed intriguing connections 
to the black-hole thermodynamics ^ ^ ^. 

Monte Carlo simulation of \wg&-N gauge theories is expected to be very useful in order to con- 
firm these conjectures or to make use of them. Indeed, the totally reduced models (the gauge 
theory reduced to D = dimension) have been studied in refs. |8|]. In the D > 1 case, some sort 
of "discretization" is needed in order to put the theory on a computer. However, lattice simula- 
tion of supersymmetric gauge theories is not straightforward. In some cases the lack of manifest 
supersymmetry just necessitates cumbersome fine-tuning, but in the worse cases the chiral and/or 
Majorana nature of fermions makes it difficult to even formulate an appropriate lattice theory. We 
propose to circumvent all these problems inherent in the lattice approach by adopting a non-lattice 
approach for one-dimensional supersymmetric gauge theories. This approach, in particular, 
enables us to put M theory on a computer using the Matrix theory [|l|], which takes the form of a Id 
UCA'^) gauge theory with 16 supercharges. 

Here we demonstrate our approach by studying the same model but in a different parameter 
region, which corresponds to the 't Hooft \wg&-N limit at finite temperature [10]. In the strong 
coupling limit the model has a dual description in terms of the A'^ DO-brane solution in type 
IIA supergravity. Our results provide highly nontrivial evidences for the conjectured duality. In 
particular, the energy (and hence the entropy) of the non-extremal black hole has been reproduced 
by solving directly the strongly coupled dynamics of the DO-brane effective theory. 



2. SUSY matrix quantum mechanics with 16 supercharges 

The model can be obtained formally by dimensionally reducing lOd = 1 super Yang-Mills 
theory to Id. The action is given by 

S=^j^dttV \^\{D,Xif - ^-[Xi,XjY + \waD,Wa - \wa{yi)afi [X;, Wfi] } , (2-1) 

where Dt = df — i[A{t), ■ ] represents the covariant derivative with the gauge field A{t) being an 
N xN Hermitian matrix. This model can be viewed as a one-dimensional U(A^) gauge theory with 
adjoint matters. The bosonic matrices X,(f) (/ = 1, • • • ,9) come from spatial components of the 
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lOd gauge field, while the fermionic matrices (a = 1, • • • , 16) come from a Majorana-Weyl 

spinor in lOd. The 16 x 16 matrices ji in ( |2.l| ) act on spinor indices and satisfy the Euclidean 
Clifford algebra /y} = 25;^. We impose periodic and anti-periodic boundary conditions on the 
bosons and fermions, respectively. The extent j3 in the Euclidean time direction then corresponds 
to the inverse temperature j3 = 1 /T . The 't Hooft coupling constant is given by A = g^N, and the 
dimensionless effective coupling constant is given by A = X/T^ . Without loss of generality we set 
A = 1, hence low (high) T corresponds to strong (weak) coupling strength, respectively. 



3. Non-lattice simulation for SUSY matrix quantum mechanics 



We fix the gauge by the static diagonal gauge 



1 



A{t) = -diag(ai,---aiv) 



(3.1) 



where can be chosen to satisfy the constraint maxa(aa) — mina(aa) < 2n using the large gauge 
transformation. We have to add to the action a term S'fp = — La<fc21n [sin |, which appears 
from the Faddeev-Popov procedure. 
We make a Fourier expansion 



A' 



•hAmn 



(3.2) 



«=— A 



r=-A' 



where (0 = ^ and A' = A — 1/2. The indices n and r take integer and half-integer values, respec- 
tively, corresponding to the imposed boundary conditions. Introducing a shorthand notation 



ki-\ hkp=n 



we can write the action ( ^ ) as 5 = 5b + Sf, where 



5b =Np 



1 £ Inco-^K^] fcxf 



' n=— A 
A' 



5f = xA^iS £ 



r=-A' 



i\ rco 



CCa-CCh 



(3.4) 



It is important that we have introduced the cutoff A after fixing the gauge non-perturbatively . This 
is possible only in Id. In higher dimensions, the momentum cutoff regularization inevitably breaks 
the gauge invariance. In the bosonic case, we have checked explicitly [^] that the results of the 
non-lattice simulation agree with the results of the lattice simulation in the continuum limit. 

Note that our action is nothing but the gauge-fixed action in the continuum except for having 
a Fourier mode cutoff. This leads to various advantages over the lattice approach proposed in ref. 



[11]. Supersymmetry, which is mildly broken by the cutoff, is shown (in Id Wess-Zumino model) 
to be restored much faster than the continuum limit is achieved. In fact, the continuum limit is also 
approached faster than one would naively expect from the number of degrees of freedom. These aie 
understandable from the fact that the modes above the cutoff are naturally suppressed by the kinetic 
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term. A further (albeit technical) advantage of our formulation is that the Fourier acceleration, 



which eliminates the critical slowing down completely 012[], can be implemented without extra 
cost since we are dealing with Fourier modes directly. We consider that all these merits of the 
present approach compensate the superficial increase in the computational effort by the factor of 



0(A) compared to the lattice approach [ |11[ ] with the same number of degrees of freedom. 

The fermionic action may be written in the form Sf = ^-y^Aar-BpsVarVps' where we have 
expanded \j/ar = La=i Var^^ terms of U(A^) generators t^. Integrating out the fermions, we obtain 
the Pfaffian Pf^, which is complex for generic configurations of the remaining bosonic variables. 
However, it turns out to be real positive with high accuracy in the temperature region studied in the 
present work. Hence we can replace it by |Pf^| = det(^^/'^), where & = One can then 

apply the Rational Hybrid Monte Carlo algorithm [ [l3| ] to study the system in an efficient way. 

4. Results 




Figure 1: The energy (normalized by A^^) is plotted against T. The dashed line represents the result obtained 
by HTE up to the next leading order for N ~ 12 Jfj]]. The solid line represents the asymptotic power-law 
behavior at small T predicted by the gauge/gravity duality. The upper left panel zooms up the region, where 
the power-law behavior sets in. 

In fig. ^ we plot the internal energy defined by E = ^(j3^), where ^ is the free energy 
of the system. Our results interpolate nicely the weak coupling behavior — calculated by the high 
temperature expansion (HTE) up to the next leading order [|l^] — and the strong coupling behavior 
= 7.4 • predicted by the gauge/gravity duality from the dual black-hole geometry Jl^]. 
The power-law behavior sets in at T ~ 0.5, which is reasonable since the effective coupling constant 
is given by A = 1/r^ in our convention. 

In ref. [Q] the Gaussian expansion method was applied to the present model, and the energy 
obtained at the leading order was fitted nicely to the power law E /N^ = 3.4 • T^-' within 0.25 < 
1. This result is in reasonable agreement with our data at T ~ 1, but disagrees at lower temperature. 

In fig. ^ we plot the absolute value of the Polyakov line P = ^Y^a^i^'"" , which is the order 
parameter for the SSB of the U(l) symmetry. It changes smoothly for the range of T investigated. 
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Figure 2: The Polyakov line is plotted against T . The dashed line represents the result of HTE up to the 
next leading order for =12 [|l^]. The dotted line represents a fit to eq. (4-.1 ) with a = 0.15 and b = 0.072. 




Figure 3: The Polyakov line for the bosonic model |]16[]. The dashed line represents the result of HTE up to 
the next leading order for N ~ 16 



which implies the absence of a phase transition as predicted by the gauge/gravity duality ||]. At 
low T it can be fitted nicely to the asymptotic behavior characteristic to a deconfined theory: 



{\P\)=exp 



+ 1 



(4.1) 



This is in striking contrast to the bosonic case [16] shown in fig. ||for comparison^. 



'in ref. [ |lq ] it was found that there are actually three phases in the bosonic model. The intermediate phase appears 
in a very narrow range of temperature Tq2 <T <Tqi, where T^i = 0.905(2) and 7^2 = 0.8761(3), and it is characterized 
by the non-uniform eigenvalue distribution of the holonomy matrix. The order of phase transitions are second order at 
r = Tci , and third order at T = 7^2- 
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5. Summary and future prospects 

We have presented the first Monte Carlo results for the maximally supersymmetric matrix 
quantum mechanics. The non-lattice simulation enabled us to study the low temperature behavior, 
which was not accessible by HTE. This provided highly non-trivial evidences for the gauge/gravity 
duality. In particular, we observed that the internal energy asymptotes nicely at low temperature to 
the result obtained from the dual black-hole geometry. 

Our results suggest that not only the power but also the coefficient of the power-law behavior is 
reproduced correctly by the gauge theory in the N and X —^oo limits. This implies that we were 
able to identify the microscopic degrees of freedom, which accounts for the Bekenstein-Hawking 
entropy for the lOd non-extremal black hole. They are nothing but the open strings attached to the 
DO-branes, which are described by the gauge theory. This should be compared with ref. [|^], which 
studied extremal black holes and relied on the supersymmetric non-renormalization theorem. 

Assuming the duality to hold in the stronger sense, one may go on and investigate the quantum 
and stringy corrections to the black-hole thermodynamics from the gauge theory side as finite-A'^ 
and finite-A effects. In particular it would be interesting to understand the physical meaning of the 



infrared instability observed in our simulation [10| from that perspective. 



When we simulate M theory, we should impose periodic boundary conditions on fermions, 
and then the system with finite N corresponds to a sector of M theory compactified on a light-like 



circle [18|. However, the Pfaffian will not be close to real positive due to the fermionic zero modes 
unlike the situation in the present work. That may cause a technical problem known as the sign 
problem when one tries to investigate the \aige-N behavior. The same problem occurs in Monte 



Carlo studies of the totally reduced models, in which the phase of the Pfaffian is speculated 019| ] 



to induce the spontaneous breaking of SO(IO) symmetry down to S0(4). This pattern of SSB 



is indeed suggested by the Gaussian expansion method |2C], and it provides a natural scenario 



for the dynamical generation of 4d space-time In ref. [g^] it was conjectured, based on the 



Eguchi-Kawai equivalence, that a similar phenomenon occurs also in the supersymmetric matrix 
quantum mechanics. We hope to address such an issue from first principles by using the non-lattice 
simulation method together with the idea proposed in ref. ^ to overcome the sign problem. 
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